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Abstract 

An asymmetric exclusion model on an open chain with random 
rates for hopping particles, where overtaking is also possible, is stud- 
ied numerically and by computer simulation. The phase structure of 
the model and the density profiles near the high and low density coex- 
istence line are obtained. The effect of a site impurity is also studied. 
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1 Introduction 



Over the past few years, a great deal of progress has been achieved in our un- 
derstanding of boundary and disordered-induced phase transitions in driven 
diffusive systems. Our understanding of such macroscopic systems, which 
are far from equilibrium, is still in its infancy since we need in this case, to 
study in detail the stochastic macroscopic dynamics of the system. It is only 
in one dimension that various exact and approximate solutions (See |l|, || for 
a review) have helped us to study the macroscopic stochastic dynamics of 
model systems, i.e. by writing the master equations of the underlying Markov 
process in suitable coordinates and solving these equations either exactly or 
by numerical-simulation techniques. The simplest such model is the Totally 
Asymmetric Simple Exclusion Process (TASEP), first introduced by Spitzer 
in 1970 0, as an example of an interacting stochastic process [|, |], |(|. This 
is a model incorporating a collection of random walkers, hopping with equal 
rates to their right empty site on a one-dimensional lattice and interacting 
via exclusion. There is now a rather extensive literature on ASEP and its 
various generalizations (see ]T|, |2| and references therein). 

Among the generalizations, we mention one which is relevant to the 
present paper, namely the introduction of disorder in the hopping rates of 
particles J?J [LC], [11], |13| . In []7], §, a closed system is considered, there is 
no injection to or extraction of particles out of the system. Each particle p , 
has a hopping rate p^, taken from a random distribution, the particles keep 
their initial order and no overtaking is possible. The bulk density of particles 
p, is an important control parameter, and it is shown that for p > p c , where 
p c is a critical density, a segregated phase appears behind the slowest particle 
on the ring, a phenomenon shown |7[] to be similar to Bose condensation in 
equilibrium statistical physics. 

The same process has been simulated |]T3| on an open chain where each 



particle enters the left end of the system with rate a and leaves the right end 
with rate (3, and the intrinsic hopping probability p^ G [c, 1] of a particle p, 
is drawn from the distribution 

/(p) = tt^tt ■ (p - c r (!) 

Using mostly the values n = 1 and c = |, it has been shown in Jl3[ that the 



generic form of the phase diagram of the single species (pure) ASEP, prevails 
also in this case; that is depending on the value of alpha and beta, three 
phases develop in the system, namely low-density, high-density and maximal- 
current phases. Compared to the pure case, there is only a slight difference in 
the shape of the coexistence line between the low and high density phases. 



2 



The small effect of disorder on the phase diagram of the pure ASEP has been 
attributed [13| to insufficient time for density inhomogeneities (platoons) to 
develop up to the scale of system size. This in turn has been attributed to the 
average passing time of particles (r ~ system size L) and the approximate 

n+1 

rate of growth of platoon sizes ~ which is small compared with 

L, when L — > oo, for any n. 

The above results all apply to a model in which particles can not overtake 
each other. The appearance of high density phase and traffic jams are quite 
expected in such a model. One can however introduce models in which 
particles can overtake each other ||. It's natural to expect that fast particles 
have a possibility to overtake slow ones, and one may ask to what extend this 
new element may affect the nature of steady state and the phase structure of 
the system. The model of admits an exact solution |Q| by way of matrix 
product ansatz , where the phases can be determined, although the density 
profiles are still very difficult, if not impossible, to obtain analytically in this 
model. Also it is difficult to obtain any analytical result concerning other 
effects like the effect of impurities or blockage. 

In this paper we study the model of || by computer simulations, and will 
obtain the phase structure of the model and density profile near the high- 
density /low- density coexistence line, and compare our results with those of 
flTB| . We also study the effect of a fixed blockage in a simple two-species 
model. 



2 The model 

The model which has been first introduced in || and shown to be exactly 
solvable via the matrix product ansatz, refers to a process in continuous time, 
in which each particle of type fi, has an intrinsic hopping rate to its right 
empty site, as in the models considered in |7|, |8|, [OJ, but when this particle 
encounters a site already occupied by a particle of type fjf (with ly < v^) , the 
two particles exchange their sites with rate — ty , as if the faster particle 
stochastically overtakes the slower one (Figure 1). 

The need for exact solvability of the model, enforces the injection and 
extraction rates of particle \x to be fixed as 

«m = —=r (2) 



and 



pv 



P, = ^+(3-l (3) 
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Figure 1: The elementary processes. 

where p is the number of type of particles and 

1 P 

^:=~yV ( 4 ) 
p /-^ 

It is assumed that there are p type of particles with hopping rates v\ < v 2 < 
V3<---<v p . 

Here a and (3 are representing the total injection and average extraction 
rates and j3 > 1 — where v\ is the slowest hopping rate. Note that for 
each particle fj, , is its velocity in an otherwise empty lattice, averaged 
over many realizations of the process. Thus, hereafter, we use the words 
intrinsic hopping rate and intrinsic velocity interchangeably. The model can 
be formulated for continuous distribution of intrinsic hopping rates which 
we denote by o~(v). In this case the rate of injection of particles of intrinsic 
velocity between v and v + dv is given by aa(v)~dv, where the extraction 



rate is given by (= + (3 — l)dv. By exactly solving this model [|LT], using the 
matrix product ansatz, the phase diagram of this model was determined. It 
was shown that besides a and (3, the phase diagram depends crucially on a 
characteristic of the distribution a(v), denoted by l[a] and given by 

t> 2 vv 
v{ [v — V\Y 

where v\ is the slowest velocity and the average is taken with respect to 
the distribution a{y). 

If a(v) is such that l[o~] < 0, then the phase diagram is similar to the pure 
ASEP, although the shape of coexistence line and the position of the triple- 
point do have essential dependence on a. However when l[a] > 0, then there 
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is no high-density phase in the diagram (Figures 2a and 2b). Specifically for 
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Figure 2: Typical phase diagrams for two values of /. 

LD low density HD high density MC maximal current 



distribution of the type(l), the size and in fact the existence or non-existence 
of the high-density phase depends crucially on the parameters c and nfl. This 
means that if the chance of entrance of slow particles into the system (cars 
into a highway) is sufficiently low, and if there is a possibility for overtaking 
slow cars, there will not be a high-density phase (a traffic jam), regardless of 
the value of a and (3. Note that in any case (3 is restricted by (3 > 1 — |, 
so that the extraction rates of all types of particles are positive. In case 
1 — '= < (3 < 1 — the steady state will be trivial and the system will 
eventually be filled with particles of type 1, since in this case all types of 
particles except type 1 are extracted from the system. In any other case (e.g. 
1 — ^ < /3 < 1 — =■), the steady state depends on the initial conditions, since 
neither particles of type 1, nor type 2, are extracted and the steady state 
depends on the order of injection of these particles. 

If to each site k, we assign a random variable r M (fc) which takes the value 
1 only if this site is occupied by a particle of type /i, then the average density 
of particle of type p at site k, denoted by p M (fc) is 



(6) 



The exact criteria is given in the main text. 
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where ( ) means the long time average, in the steady state. Note that in 
accordance with [13] and in contrast to the periodic lattice |7], ||, a separate 
disorder average is not necessary, since new particles are constantly injected 
into the system. Two other quantities are of interest, defined as: 



m ._ 



1 L 



*(*)> 



(7) 



k=l 



which are respectively the density of particles of type p, average over the 
lattice sites; and the total local density of particles, irrespective of their 
type. One can also define the current of type fx of particles, which in the 
steady state is independent of site and in view of the definition of the process 
is given by || 

J" = vM_ x e k ) + - v)«.i^> " £(*V - v ^ti<) (9) 

where e*, := 1 — X/u=i T k ■ Note that (e^) stands for the probability of site k 
being empty. 

The independence of this current of the site k, and its equality to the 
input and output current of particles of type fi, given by 

JL = ^(ei) (10) 



■C = (0 + f -i)<r£> (ii) 

gives us a set of equations, which in the mean field approximation 
(T^(k)T u (k + 1)) ~ (T^(k))(r u (k + 1)) enables all the local densities ff(k) 
and hence all the currents J M , to be evaluated numerically. 

We also study the effect of a site impurity at a given site, say k, by 
multiplying all the hopping rates by A (A < 1), when computing the currents 
of particles on the links (k — l,k) and (k, k+1). The densities and the current 
can again be computed for this case, for various A's and Fs. 
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3 Numerical solutions and simulation results 



We consider a lattice of 200 lattice sites and simulate the following process 
with random updating on this lattice: 

a) A two species process, where the hopping rates are taken to have the 

distinct values v% and v 2 , with equal probabilities. 

b) a multi-species process, where the hopping rates are taken from the dis- 

tribution 

^) = 7 I ^TT(--^ > o<v<l (12) 



c) a two species process, where there is a fixed impurity on one of the sites 
of the lattice. When the particles reach this site, they slow down their 
hopping rate by a factor A. The position of the impurity is denoted by 
k, where < k < L. 

In cases a and c, we complement our simulations by a mean field solution 
obtained by numerically solving the equations «/£ = ••• = = ■ ■ ■ = J^ ut . 



Results 



a) A two species process 

Fixing the values of v\ and v 2 , we obtain the bulk density p, for different 
values of (3, as function of the injection rate a. 

The results are shown in Figures 3a and 3b, for the two different choices 
of hopping rates. 

It is seen that for every (3 , there is a jump discontinuity in the bulk 
density at a critical value of a c = a c (j3). The value of the jump vanishes at 
a certain point (a c ,/3 c ), which makes the triple point in the phase diagram. 
The corresponding phase diagrams are shown in Figure 4. 



b) A multi species process 



In this case we consider a continuous distribution of hopping rates. For the 
sake of comparison with fl3fl , we take the distribution (12). 



For this distribution l[a] is calculated to be JTO 
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/ -oo if 0<n<l 

[al ~ I (I) 2 - (n + + (^TXT^) if 1 < n ^ 

where v := f c va(v)dv = c + — c). 

For n = 1, we expect from the above formula that, since l[a] < 0, we have 
the usual three phases seen in the pure case, namely the low density, the high 
density and the maximal current phases. this expectation is borne out in our 
simulation, where we have obtained the average density as a function of a, 
for various values of ft (Figure 5a). 

As far as c = ~, one can see after a simple calculation that l[a] can never 
be positive. To obtain a positive we take n = 2 and c = | (a broader 
spectrum of hopping rates) for which / is found to be equal to 3. In this case 
the results of simulations are shown in Figure 5b. 

It is seen that when l[a] > 0,no jump discontinuity exists in density, 
implying that there is no high density phase in the phase diagram. For all 
ft, there is a critical a c independent of ft, after which we enter the maximal 
current phase ; this is in accord with the exact a analysis of [IT|. The 
corresponding phase diagram is shown in Figures 6a and 6b. 



c) The effect of a site impurity 

In this problem, our aim is to see what effect a site impurity has on the density 
profiles of particles, when the particles have the possibility of overtaking each 
other. For simplicity, we consider a two species model. 

It is expected that such a blockage will mostly affect slow particles, since 
fast particles have a chance of passing by the particle accumulated behind 
the blockage. This is indeed seen in our numerical solution. Figure 7a and 
7b depict the density profiles of both types of particles (fast and slow) for 
two different values of V\ and v<i . Density profiles changes versus the strength 
of impurity A (the slowing down factor) is shown in Figure 10. 
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(b) vi = 0.8, v 2 = 1.2 

Figure 3: The bulk density as function of a for various values of j3, in the 
two species model. 
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Figure 4: The phase diagram for problem a. 
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(a) n = 2, c = | 
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jure 5: Bulk density versus ct for the distribution (12)(type number 
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jure 6: The phase diagram for problem b (type number = 200). 
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Figure 7: Density profile of slow and fast particles in the presence of an 
impurity (A = 0.5, a = 0.8, (3 = 0.8). 
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Xx 10 

Figure 8: Total density profile versus the strength of impurity (v\ = 0.8, Vi = 
1.2, a = 0.8,(3 = 0.8). 
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4 Conclusion 



In this paper we have studied an asymmetric exclusion model on an open 
chain with disorder in particle hopping rates. Following the method of 



13], and guided by the analysis of Wu PL we have shown how the distribu- 



tion of hopping rates affects the generic phase diagram of the single species 
ASEP. In particular we have taken the hopping rates from the distribution 



<v)= ^l^ iv-cY^cKvKl (14) 



also studied in ||13|| , and have shown that for n = 2, the high density phase 
may or may not exist depending on the value of c; i.e. the width of the 
distribution. We have also studied how in a single two species process, with 
hopping rates V\ and V2, a single site impurity affects the density profiles of 
fast and slow particles. 
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